Oceanic coastline and super-universality of percolation clusters. 
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New fractal subset of a rough surface, the "oceanic coastline", is defined. For random Gaussian 
surfaces with negative Hurst exponent H < 0, "oceanic coastlines" are mapped to the percolation 
clusters of the (correlated) percolation problem. In the case of rough self-afnne surfaces (H > 0), 
the fractal dimension of the "oceanic coastline" d c is calculated numerically as a function of the 
roughness exponent H (using a novel technique of minimizing finite-size effects). For H — 0, the 
result d c w 1.896 coincides with the analytic value for the percolation problem (91/48), suggesting 
a super-universality of d c for correlated percolation problem. 
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Random surfaces are common objects which can be 
found in various situations; good examples are deposited 
metal films Q , fractured surfaces 0, U , streamfunctions 
of two-dimensional (2D) turbulent flows Q etc. For var- 
ious applications, it is important to know the statis- 
tical geometrical properties (statistical topography) of 
such surfaces. This problem has been studied with a 
considerable success, both numerically and analytically 

Here we introduce a new scale-invariant measure of 
such surfaces, the fractal dimension of the "oceanic coast- 
lines" . An "oceanic coastline" includes the outer perime- 
ter of the ocean (the "continental coastline"), and the 
coastline of "oceanic islands" , i.e. of such islands, which 
are adjacent to the largest (possibly infinite) connected 
flooded region, when water floods the surface up to a 
certain level; "oceanic islands" are to be opposed to such 
islands, which are inside the lakes. Thus, each flooding 
level defines its own "oceanic coastline". Here we show 
that for the so called "monoscale" random surfaces (with 
bounded height, cf. [5(), the "oceanic coastlines" can be 
mapped to the percolation clusters of the correlated per- 
colation problem. For rough surfaces (with unbounded 
height), the fractal dimension d c is calculated numeri- 
cally, using a novel technique of minimizing the finite-size 
effects. The results suggest that for the correlated per- 
colation (and hence for the "monoscale" surfaces), there 
is a super-universality, d c = 91/48. 

An important characteristic of random surfaces [given 
by the surface height ip[x,y) = ip(r) over the two- 
dimensional plane] is the Hurst exponent H; in the case of 
Gaussian surfaces, all the statistical topography scaling 
exponents are functions of H. The most intuitive defini- 
tion of the Hurst exponent is given by the scaling law of 
the surface height drop at a distance a = \a\ (cf. 

(^(r)-^(r + a)) 2 )oc\a\ 2H . (1) 

Here angular braces denote averaging over an ensemble 
(or over the radius vector r); we assume that a is larger 
than the lower cut-off scale ao, and < H < 1. More 
general definition (not bounded to < H < 1) of the 
exponent H is given via the power spectrum P& [with 
(ipkipk') = Pk8(k + k'), where -0fc is the Fourier transform 



of the surface height; cf. Q]: 

P k (x\k\- 2H - 2 for Ifcl^a^ 1 . (2) 

According to the value of H, surfaces can be divided 
into three classes 

(i) H < 0: "monoscale", surface looks like a landscape 
with hills of approximately same size, the surface height 
is bounded; 

(ii) < H < 1: "rough" , looks like rocky mountains with 
smaller and higher peaks and unbounded height; 

(Hi) H > 1: "smooth", the largest scale defines the slope 
of the surface. 

For "monoscale" surfaces, there is a simple mapping 
between statistical topography andpercolation problem, 
developed by Ziman and Weinrib [UGH [the minima of 
the surface define the sites of an irregular percolation 
lattice; the saddle-points S%j (between i-th and j-th min- 
ima) correspond to the bonds of the lattice (which can be 
drawn as the steepest descend paths from Sij to the i-th 
and j-th minima); when the surface is flooded by water, 
a bond is declared broken if the corresponding saddle is 
("dry"), i.e. above the water level; the water level itself is 
the counterpart of the bond- keeping probability p] . This 
mapping can be also applied to the case of rough sur- 
faces. Then, however, the lattice problem is no more the 
(correlated) percolation problem. In particular, there is 
no critical probability p — p c at which the global connec- 
tivity disappears. Indeed if we "look far enough" , we can 
find as deep (hence also as wide) connected wet region 
as we want, because the height is unbounded. Therefore, 
the mapping between the percolation problem and statis- 
tical topography allows us to extend various concepts of 
the percolation problem outside the range of applicability 
of the very percolation problem. 

A straightforward example here is the hull of a perco- 
lation cluster: this is defined as the outer perimeter of a 
percolation cluster. Since the percolation cluster corre- 
sponds to a connected wet region, its outer perimeter cor- 
responds to an isoline of the surface. It has been shown 
that for uncorrelated percolation, its fractal dimension 
dh = 7/4 [Tl|| . Note that percolation with weak correla- 
tion (with H < —3/4) belongs to the same universality 
class as the uncorrelated percolation Meanwhile, for 
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rough surfaces, dh ~ 1.5 — 0.51? 0, with an exact result 
d h = 1.5 for H = @. Therefore, for -3/4 < H < 0, 
one can expect a non-trivial decreasing function dh(H), 
with d ft (-3/4) = 1.75 and c4(0) = 1.5. Further, for the 
percolation problem, one can distinguish between the hull 
and the unscreened perimeter. This is because the hull 
is very twisted, and typically, there are narrow "necks" 
of peninsulas separating "gulfs" . We can postulate that 
these "necks" are too narrow, so that the coastline of the 
respective peninsulas is "screened" ; the rest of the coast- 
line is called "unscreened perimeter" . For percolation 
clusters, the fractal dimension of the unscreened perime- 
ter is d u = 4/3 Since for rough surfaces, narrow 
"necks" are untypical, d u {H) = dh{H) for H > 0. Com- 
paring d u {— 3/4) = 4/3 and d u (0) — 1-5, one can expect 
an increasing function d u (H) for — 3/4 < H < 0. 

Now we are ready to turn to the percolation cluster. 
Each bond of the cluster is touched by the perimeter 
(hull) from both sides; thus, the length of the net perime- 
ter scales as the cluster size. However, the perimeter 
touching it can be either external or internal, see Fig. 1. 
The sum of external and internal perimeters is the sum 
of (a) the coastline between the "ocean" (=cluster) and 
the "continent" (=infinite cluster of broken bonds), and 
(b) the coastline between the "ocean" and the "islands" 
(=broken bond clusters inside the cluster of undamaged 
bonds); we call this the "oceanic coastline". Note that 
the coastlines of the "lakes" inside the "islands" , "is- 
lands" inside the "lakes" inside the "islands", etc., do not 
touch the "ocean" , i.e. the bonds of the cluster, and are 
not relevant here (including these coastlines, too, would 
result in the set of all the isolines at the given height, the 
fractal dimension of which df = 2 — H for < H < 1; 
for H < 0, there is a super-universality, df =2 ). 

In order to calculate numerically the fractal dimension 
of the "oceanic coastlines" , we have opted for the follow- 
ing procedure. First, 1+1-dimensional [(1+1 )D] random 
surfaces (which are assumed to belong to the same uni- 
versality class as the statistically isotropic 2D surfaces) 
were generated on the lattice of the four-vertex model 
|13| . The one-dimensional discretized fractional Brown- 
ian functions Jh (i) were generated as fractional integrals 
of uncorrelated random sequence of "spins" Si = ±1 as 
follows. The aim function gn(i) is found as 



J=0 



(3) 



Further, the discretized function (*) is found incremen- 
tally, starting from H = 0, by rounding <?#(«) to such a 
nearest integer that 



Mi) - Mi - 1) = ±1- 



(4) 



Finally, the (1+1)D function is given by the superposition 
of two ID-functions. 




FIG. 1: All the bonds of a percolation cluster (black bold 
lines) are touched from both sides, either by external perime- 
ter (hull, dark grey lines) or internal perimeter (light grey 
lines). 



After generating a 2D surface on a square polygon of 
side length N, we started "flooding" the polygon. To this 
end, at the perimeter of the polygon, the lowest point 
Pq was found. Further, "water" was incrementally "in- 
jected" into that point, until the flooded region connected 
opposite edges of the polygon (either left and right or 
top and bottom), see Fig. 2. Such a procedure guaran- 
tees that only the "ocean" (and not the "lakes" ) becomes 
wet. At that stage, three quantities were recorded: L\ — 
the coastline length, — the number of cells touching 
the coastline, and L3 — the coastline length immediately 
before achieving the critical flood level. The procedure 
was repeated (for each polygon size and Hurst exponent 
value) 10 s or more times. The covered polygon sizes are 
given by formulae N = 128 • 2 k and N=196 • 2 fe , with 
k = 1, 2, 3,4, and N — 4096 (the largest polygon size was 
used only for H = 0). The simulations were performed 
on the cluster of ten AthlonXP-1700 workstations during 
three months, using an assembly-optimized code. 

It was expected that the finite-size scaling of the aver- 
age values k(N) = (Li(N)) allows us to estimate the frac- 
tal dimension d c . Such an asymptotic behavior was ob- 
served, indeed, for larger Hurst exponents (H > 0.375). 
However, for smaller values of H (and in particular, for 
H = 0), the convergence ^(A^) — > AN dc was very slow, 
see Fig. 3 (here A denotes a constant). Consequently, if 
we assume that U can be expanded asymptotically, 



F H {i,j) = M(j)-fm(j). 



(5) 
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the non-leading exponents (y ^ 0) have to be very 
close to the leading term. In order to resolve this diffi- 
culty, we made an assumption that for the three recorded 
quantities U (i — 1,2,3), the three leading exponents 
{y = 0, 1, 2) are the same, a U i — a^i. The validity of this 
assumption is evident in the case of the leading expo- 
nent a„o = cZ c , but not so clear for the other exponents. 
However, if we accept the assumption, it can be verified 
later using the numerical results. Thus, if we assume that 
a vi = O-jiii one can construct such a linear combination 
of the three quantities Zj, for which two terms (out of 
the three leading ones) cancel out. It suffices to find a 
root-mean-square fit by minimizing the expression 



S(d) = 



N d (i)-A 1 l 1 (i)-A2l2(i)-A 3 h(i) 



Here, Oi is the variance of the expression A\l\{i) — 
^-2h(i) ~ and can be expressed via the indi- 

vidual variances of the quantities Z(i), and the pair- wise 
correlation coefficients; all this is easily calculated, using 
the Monte-Carlo simulation results. The resulting curves 
\og 10 [S (d) / (n — 4)] are given in Fig. 4; here n is the num- 
ber of terms in the sum S(d) (so that n — 4 is the number 
of degrees of freedom). Note that the two smallest data- 
points (with Nq = 128 and N\ — 196) have been excluded 
from the sum, otherwise, the residual squares would have 
been unacceptably large [resulting in S/(n — 4) > 1]. 
Evidently, for N < 196, the higher (y > 3) terms of 




FIG. 2: The polygon (N = 2048, H = 0.5) has been flooded 
so that the flooded region (grey area) connects left and right 
edges of the polygon. The striped patterns (in horizontal and 
vertical directions) are due to the (1+1)D geometry of the 
surface. 
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FIG. 3: Finite-size effects for H = give rise to a very 
slow convergence of the estimated fractal dimension d c (N) = 
ki[h(N k+1 /h{N k )]/ki(N k+1 /N k ), N = y/jV k N k+1 . 



the expansion (6) become important. The existence of 
three sharp minima for all the curves in Fig. 4 confirms 
the assumption which was made about the equality of 
the leading exponents. Indeed, if the second and third 
asymptotic expansion terms of Zi, 1%, and Z3 were linearly 
independent, they would not cancel out in Eq. (7), and 
there would be no minima with S/ (n — 4) ~ 1. 

This novel method of interpreting the finite-size simu- 
lation results allowed us to calculate not only the fractal 
dimension of the "oceanic coastlines" (see Fig. 5), but 
also the second and third exponents in the asymptotic 
expansion of the length of "oceanic coastlines" (Fig. 6). 




FIG. 4: The quantity log 10 [S (d) / (n - 4)] [see Eq. (7)] is plot- 
ted versus d for different values of the Hurst exponent H 
(each curve is labeled with a number indicating the value of 
H). For all the curves, three minima are clearly distinct. At 
these minima, d = on, where Oi is the i-th leading exponent 
of the asymptotic expansion of the quantities Zi, h, and Z3. 
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The values d c = 1.8965 ± 0.0025 and a x = 1.58 ± 0.01 
for H = are relatively close to each other, exactly as it 
was expected, judging by the finite-size scaling in Fig. 3. 

Note that for H < —3/4 (uncorrelated percolation), 
there is analytic result d c = 91/48 w 1.8958 which 
is indistinguishable (at our numerical uncertainty) from 
our result for H = 0. While the author is not able to 
give a rigorous proof, it seems natural to assume that 
the function df is a monotonically decreasing function of 
H . Indeed, a larger H means a stronger influence of long 
scales, which results in less twisted isolines, and a lesser 
number of small islands (note that these arguments do 
not apply for d u : excessive twists of the hull can effi- 
ciently screen and smooth the perimeter). This leads us 
to the conjecture d c = 91/48 for H < 0. 

In conclusion, we have introduced the concept of 
"oceanic coastlines" , and shown that they are the exten- 
sion of the concept of percolation clusters. The finding 
that their fractal dimension d c is clearly distinct from 
both the dimension of single isolines d^ and the dimen- 
sion of the full set of isolines at a given height dj is partic- 



4 



\ 0.2 



0.4 



0.6 



0.8 



H 



0.01 

[-0.02 



t/ 



/ 



/ 



/ 



-1 A 
i ,*+ 




0.4 


0.6 


T T 

0.8 H 


1,2 

-1 


0.2 X 








0,8 










0,6 




^2 






0,4 










0,2 











FIG. 6: Second and third largest exponents (qi and 02) of 
the asymptotic expansion of the "oceanic coastline" length 
are plotted versus the Hurst exponent H. 

ularly important, because in the context of Earth's land- 
scapes, the coastlines of "oceanic islands" have been often 
mistakenly treated as the full set of isolines @. Based 
on the numerical results, we have conjectured that for 
H < 0, the fractal dimension of percolation clusters is 
super-universal, d c = 91/48. Last but not least, a novel 
and universal approach to the problem of interpreting the 
finite-size simulation results has been devised. 



FIG. 5: The mismatch between a linear estimate du n = 21 _ 
H|| and the fractal dimension of the "oceanic coastline" d c 
is plotted versus the Hurst exponent H . 
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